Soliton phase near antiferromagnetic quantum critical point in Q1D conductors 
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In frameworks of a nesting model for Q1D organic conductor at the antiferromagnetic (SDW) 
quantum critical point the first-order transition separates metallic state from the soliton phase 
having the periodic domain structure. The low temperature phase diagram also displays the 2nd- 
order transition line between the soliton and the uniformly gapped SDW phases. The results agree 
with the phase diagram of (TMTSF)2PF6 near critical pressure [T. Vuletic et al., Eur. Phys. 
J. B 25, 319 (2002)]. Detection of the 2nd-order transition line is discussed. We comment on 
superconductivity at low temperature. 
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We address the issue of quantum critical point (QCP) 
for the itinerant antiferromagnetism (AFM) (also known 
as spin-density wave state (SDW)) in Q1D compounds. 
As it was pointed out in |2j , the standard renormalization 
group (RG) analysis is not applicable for superconduc- 
tivity QCP and theoretical models for AFM with nesting 
features because no expansion in a form of a Landau func- 
tional is possible near QCP at T = in these cases. QCP 
with increase in pressure have been observed for SDW in 
the Bechgaard salts (TMTSF) 2 X (with X=PF 6 [Ik AsF 6 
0) and recently investigated in many details inQ. Ac- 
cording to [jj], the phase diagram of (TMTSF) 2 PF 6 near 
its critical pressure, p c ~ 9.4kbar, indeed, has a rather 
complicated character. However, the experimental fact 
that the transition from metallic to SDW state is of the 
first order suggests that quantum fluctuations do not play 
decisive role near p c , and the overall phase diagram near 
p c for (TMTSF)2PF6 can be understood already on a 
mean field level. 

The most unexpected finding in [l| is the discovery of a 
pressure interval below p c , where SDW and metal phases 
(SDW and superconductivity (SC) at lower T) coexist as 
parallel domains running perpendicular to the chain di- 
rection. Such coexistence is difficult to understand for the 
1st order transition that takes place at constant pressure 
0. In what follows we propose that such domains can 
be interpreted as formation of a new phase of the soliton 
walls suggested first theoretically for the Q1D electron 
spectrum in 0, 0, ■ Domains inside the SDW interval 
were also observed in ||. We shall see that it is rather 
natural to expect SC in these domains at lower temper- 
ature. In this presentation we restrict ourselves mainly 
by the effects related to the physics of SDW. Supercon- 
ductivity appears at lower temperature than the onset of 
the SDW/metal phase transition, and we only comment 
on it at the end. 

For the SDW description we adopt the standard Q1D 
model of two open Fermi surface (FS) sheets with the 
energy spectrum of free electrons in the form: 

e(p) = ±v F (p x TPf) +t±(jP±), (1) 
with dispersion in the direction perpendicular to the 



chains 



t±(p±) =t b (p±)+t' b (p ± ). 



The condition 



(tb(Px)) Vl =(t' b (Px)) Vl =0 



(2) 



(3) 



is assumed to conserve the number of electrons. The 
angular brackets mean averaging over p± : 



< .. > P± = J (..){b/2n)dp ± . 

The component tb(p±) possesses the so-called "nesting" 
feature (at t' b = 0) favoring formation of the SDW gap 
along the two sheets of the Fermi surface. For that one 
needs 



e{p + Q) = -e(p). 



(4) 



" Antinesting" term, t'Apx), destroys the perfect nesting 
condition Q . Increase of t' b simulates increase in applied 
pressure. At large enough t' b SDW disappears. 

The model assumes a bare e-e interaction of the form 



u(Q)<t i( t 2 , 



(5) 



where u(Q) has a maximum at some Qo- Interaction (JSJ 
is renormalized after summing up the ladder diagrams: 



U(Q) = 



u(Q) 



l + u(Q)II(Q,wo)' 
The polarization operator 



(6) 



/j2 
— ^G(p,w n )G(p-Q,Lj n -u) (7) 



is proportional to the staggered susceptibility. Its value is 
large near a "nesting" vector Q . Most often this vector 
is chosen as 



Qo = (±2PF,ir/b,ir/c), 



(8) 



although the majority of the results below do not change, 
at least qualitatively, at different optimal " nesting" wave 



2 



vectors Qo- For brevity, in (J2J we shall leave out the 
dependence on the third direction along the c-axis Q. 

When the denominator in © becomes zero, the 
metal phase is absolutely unstable with respect to a 
SDW(CDW) nucleation. The instability line T c (t' b ,Q) 
is given by a well-known equation: 



In 



' cO 



= Re 



iW(Q, P± ) 



- * 



(9) 

where T c0 is the temperature of the SDW onset at t' b = 
and W{Q,p±) = Tiv F (Q x - 2k F ) + t ± (p y ) + t ± (p y - Q y ). 
Eg. h as been analyzed by many authors (see, e.g. 
EE lid . Il l L H ^ l. At zero temperature Eq. © writes 
down as 



(ln(TTT c0 /jW(Q,p ± ))) =0, 



(10) 



where 7 1.781 is the Euler constant. 

With increasing pressure the value of i£ increases, and 
at some critical pressure p c the SDW ordering may dis- 
appear even at zero temperature. This point could be 
considered as QCP. However, one sees that the applica- 
tion of renormalization group (RG) analysis [2| to QCP 
in Q1D case would meet difficulties because an analytic 
expansion of Eq. (|10J) at small q = Q — Qo does not exist 
due to the condition 

The integral JTOJ is convergent at any finite Q, and 
one can find critical t' b for the general form of t±(p±). 
It does not necessarily occur at Q = Qo, but at some 
Q* = Qo + q* 0. Transition into such a new SDW 
phase with "shifted" wave vector (SDW2 phase) should 
be of the second order, and, hence contradicts to the 
observed lst-order character of the transition between 
SDW and the metal phase in 0]. 

It is known that for strictly ID models of CDW (or 
SDW) the character of excitations is different from the or- 
dinary one. Instead of electron-hole pairs with excitation 
energy rj 2A, the propagating excitations are solitons (or 
soliton " kinks" ) that cost a lower energy and come about 
together with the nonlinear reconstruction of the under- 
lying gaped ground state HI EH- (TMTSF) 2 PF 6 has 
a quarter filled band (2k F — w/2a). Unlike the Peierls 
state in polyacetylene, (CH) X , commensurability effects 
play here no essential role [lq. so that only the neutral 
spin 1/2 solitons exist Below we merely borrow the 
results from the extended literature on solitons in a ID 
chain (for a review and references, see [13 )■ Although 
this literature is devoted to the physics of CDW state, 
all the results can be transferred to the SDW case, prac- 
tically, without changes. 

The energy cost for one soliton "kink" on a single ID 
chain is 



(2/tt)A . 



(11) 



For chain packed in a crystal a tunneling, t±(p±), arise 
between the chains. Therefore, the solitons form ex- 
tended states by creating a band in the transverse direc- 
tion. Instead of independent kinks on different chains, 



one arrives at a soliton wall [5j with neutral solitons oc- 
cupying two spin states with t±(j>±) < in the band . 
The energy cost A(t±) of the wall referred to one chain 
decreases, 



(12) 



A(t ± ) = (2/tt)A - 2 / t x (p±)bdp ± /2TT, 
Jt ± <o 



and may even become negative. In other words, the spon- 
taneous formation of soliton wall may be favorable at 
large enough t±. 

Eq. Ijl2|l presents our key idea for interpreting the phase 
diagram of (TMTSF) 2 PF 6 Q: with the pressure increase 
the system first reaches a critical pressure, p c ±, at which 
A(t±) = 0, and smoothly enters into a new phase - the 
soliton phase (SP). At some higher pressure (p c ) the first- 
order transition between SP and metallic phase takes 
place. 

To describe the sequence of such transitions, one needs 
expressions for the energy of SP. The problem was solved 
in ja| for the CDW ordering and for " direct" nesting vec- 
tor Qo = (±2fc F ,0). Change from CDW to SDW needs 
no explanation. Without staying on the proof, we remark 
that the gauge transformation of the wave functions, as 
in 0, removes the large "nesting" term t b (p±) from all 
expressions below. As the result, it is only the "antinest- 
ing" term tf b (p±) that enters Eq. JEJ: t±(p±) —> t' b (p±). 

We now write down the expression for the linear energy 
density Wsp of the soliton phase [|| in the limit of large 
distance between the walls(T = 0): 



Wsi 



A 2 



2irhvp 



nA(t' b ) + nE 2 _B 1 



(13) 



where n is the walls' linear density, and the term in 
()13J) correspond to the exponentially decaying interaction 
between the walls. Indeed, in this limit of large distances 
between walls, n from 



E + /h,VF 



K 



In (4A /£;_ 



(14) 



where E + m A at n — > and K(r) is the complete 
elliptic integral of the 1st kind. From i|14[l 

E_ w 4A exp(— Ao/nTwir). 

A(t b ) is given by Eq. ifl^jl with tj_(p±) = t b (p±), and 

1 b 



B(t ± ) = 



27rA \dt' b /dpxW 



(15) 



where \dt b /dpx\o is the value of the transverse velocity at 
the four point where t b (p±) — 0. At B > crossing the 
point A(t b ) — corresponds to the second-order transi- 
tion from the "homogeneous" SDW to the soliton walls 
lattice state. Negative B < would mean an abrupt 
first-order phase transition [6j. All the results (|12I15|) 
are for T = 0. 
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We now return to Eq. As the absolute instability 
line, at higher temperature it defines the 2nd-order tran- 
sition line between the metal and SDW(CDW) phase. At 
lower temperature for most t b (p±) models it corresponds 
to some "shifted" Q* . As it was mentioned above, ex- 
perimentally pj, the metal-to-SDW transition at lower 
temperature is of the lst-order. At these temperatures 
the line defined by Eq. (0 has only the meaning of a 
"supercooling" line. 

Theoretically, the positions of critical points p c and p c \ 
depend essentially on a particular form of t' b (j)±). The 
function t' b (p±) in (TMTSF) 2 PF 6 is unknown. There- 
fore, for qualitative analysis we will apply the expressions 
(|14|l . (|15|) also in the regime of dense soliton phase. We 
compare two examples: first, the oversimplified model, 
in which t' b (j>±) has the periodic step- like shape: 

tbiP ^-\-2t b , n/b<p^<2n/b, (16) 

and the second one with the tight-binding dispersion of 
the form 



t' b ( P± ) = 2t' b cos(2 P± b) 



(17) 



We calculated the energy of the soliton phase for both 
dispersion functions. At the calculation the optimal value 
of _E_ (or n) must be found to minimize the energy l|13|) of 
the soliton phase. For the first dispersion (|16fl the domain 
wall energy A(t' b ) becomes negative at t b > t c \ w 0.32A . 
This point corresponds to the 2nd-order transition from 
"homogeneous" SDW phase into the soliton phase. Then 
we compared the soliton phase energy Wsp with the en- 
ergy of metal phase 



W n 



-(l/nhv F )(lt' b (p^f) pi 



(18) 



and found the point t c ~ 0.54Ao where these energies 
become equal to each other. This point corresponds to 
the first-order phase transition from SP to the normal- 
metal phase. We see that the interval t c — t c % w 0At c 
of the existence of SP for this particular model is rather 
large. One can show that for the model described by 
Eqs. Ijl2|l - (|15|l the function l|16(l correspond to the largest 
possible interval of SP. 

For the second dispersion relation (|17|) the similar cal- 
culation gives equal values for t c and t c \: t c = t c \ — 
A /2, i.e. in this case an interval for soliton phase is 
absent. This once more shows the high sensitivity of the 
SDW phase diagram to the particular form of t b (p±). 

In experiment [1| the difference t c — t c \ ~ 0.1i c is in- 
termediate between the cases of Eq. (fTB|) and Eq. (fTTjl . 
but closer to the second case of tight-binding dispersion 
relation (|17fl than to the step-like dispersion ]l(i|l. The 
ratio (t c — t c i)/t c rs 0.1 observed in experiment [l| can be 
easily fitted by an appropriate choice of the model func- 
tion t b (p±) (for example, by adding a fourth harmonic to 
the tight- binding dispersion (|17|) . 

In Fig. 1 we show schematically the SDW phase dia- 
gram near QCP. The slope of the second-order transition 



line near T = comes about because at T ^ the sec- 
ond term in A(tj_) decreases to account for thermal ex- 
citations inside the domain wall from the occupied part 
of band t b (p±) < 0. The temperature slope, of course, 
again depends on the ^(pjj-dependence. 

In the phase diagram of (TMTSF) 2 PF 6 has been 
studied through the resistance measurements along the 
chain direction. The lst-order character of the transition 
between metal and SDW states follows from hysteretic ef- 
fects, while the appearance of domains was derived from 
the changes in resistivity behavior near p c . We suggest 
that the line of the 2nd-order transition between the in- 
sulating " homogeneous" and soliton- wall phases could be 
detected by measuring a sharp anisotropy of conductiv- 
ity near this transition Q. Indeed, according to 0, at 
low n (large distance between the domain walls) conduc- 
tivity along the chain direction remains very low because 
of exponentially small overlap between the electron wave 
functions inside each wall. As for transverse conductivity, 
it would increase linearly with n. 

One also sees that the value of t' b itself is not the 
only parameter that may determine the dependence on t' b 
(pressure). Indeed, the expression for B(t±) i|15|) essen- 
tially depend on \dt b /dp±\o as well. The i?-term is not 
analytic at t' b — ► 0; change in its sign will immediately 
change the whole physics |6|]. 

In Fig. 1 the superconductivity, seen experimentally 
in near p c and in the metallic phase at lower tempera- 
tures (T c ~ IK ), is not shown in order not to overcrowd 
the phase diagram. Meanwhile, appearance of supercon- 
ductivity is inevitable in frameworks of nesting mean- 
field models. Indeed, the logarithmic singularity that 
enters the expression for SDW polarization operator J7J), 
would appear also for the Cooper channel |l8| • At higher 
temperature the SDW phase prevails because the bare in- 
teraction J5J was chosen to be larger for the transverse 
nesting component Q±. However, with increase of "an- 
tincsting" term t' b (pressure) the SDW transition becomes 
suppressed and SC which is not sensitive at all to any of 
the ij>-terms, finally becomes more favorable. To some 
surprise, T$c for the SC transition is experimentally not 
sensitive to the formation of the SP. We explain this fact 
qualitatively by assuming that the SC coherence length, 
£o ~ Tivf/Tsc, is considerably larger than the expected 
periodicity for the SP of the order of Uvf/^o ~ fiVp/To 
resulting in considerable Josephson coupling between the 
soliton walls. However, the problem needs some further 
analysis. 

To summarize, we have shown that in the vicinity of 
QCP in Q1D materials, such as the Bechgaard salts, one 
meets with a new soliton wall phase. As the pressure 
increases from p = 0, one first crosses the line p c ± of the 
second-order transition and enters into the periodic soli- 
ton structure with a characteristic pressure-dependent 
periodicity of the order of Tivf/Aq ~ 10 3 A. At higher 
pressure, p c , by the first-order phase transition the sys- 
tem goes over into metallic state. Details of the phase 
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FIG. 1: A schematic view of the SDW phase diagram. The solid line signify the first-order transition between the metal and 
soliton phases. The dot line stands for the second-order transition between uniformly gapped SDW phase and the SP. The 
dash line represents the line of absolute instability of metal phase toward the formation of SDW. At higher temperatures this 
line represents the second-order transition line between the SDW and metal phases. At lower temperature (dash-dot line) it 
has only the meaning of supercooling line. The lower filled bars shows the region where superconductivity appears. We do not 
calculate the superconducting critical temperature. The upper filled bar shows the high-temperature region of soliton phase 
where substantial quantitative deviations from formulas 112H - 115|I may start due to the thermal excitations in the SDW and 
soliton phases. Detailed shape of all lines is sensitive to the choice of the " antinesting" term t b (j>±) (see text). 



diagram strongly depend on material parameters, in par- 
ticular on the form of the " antinesting" tunneling in the 
electron dispersion. Coexistence of SDW and supercon- 
ductivity inside soliton walls has not been investigated. 
Qualitatively one would expect the survival of supercon- 
ductivity in the " metallic sheets" (domain walls) perpen- 
dicular to the chain direction 0. 
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